Passage of radiation through wormholes 
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We investigate numerically the process of the passage of a radiation pulse through a wormhole 
and the subsequent evolution of the wormhole that is caused by the gravitational action of this 
pulse. The initial static wormhole is modeled by the spherically symmetrical Armendariz-Picon 
solution with zero mass. The radiation pulses are modeled by spherically symmetrical shells of 
self-gravitating massless scalar fields. We demonstrate that the compact signal propagates through 
the wormhole and investigate the dynamics of the fields in this process for both cases: collapse of 
the wormhole into the black hole and for the expanding wormhole. 
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I. INTRODUCTION 

One of the most interesting features of the theory of 
general relativity is the possible existence of spacetimes 
with wormholes [E S S 0, S @|- The wormholes are 
topological tunnels which connect different asymptoti- 
cally flat regions of our Universe or different universes 
in the model of Multiverse • 

Recently the hypothesis that some known astrophysi- 
cal objects (e.g. the active nuclei of some galaxies) could 
be entrances to wormholes was considered by Kardashev, 
Novikov and Shatskiy and by Shatskiy [Iflj]. If some 
of these wormholes are traversable, then radiation and 
information, at least in principle, can pass between two 
regions of the Universe or from one universe to another 
universe in the model of the Multiverse. Furthermore, 
if a wormhole exists, in principle it is possible to trans- 
form it into a time machine (space-time with closed time- 
like curves) [111, [l^l ■ Therefore, it is very important for 
both the theory and the development of the hypothesis 
about the real existence of the wormholes to analyze the 
physics of the passage of radiation through a wormhole. 
According to the wormhole models of general relativity, 
in the absence of any special matter, the throat of the 
wormhole pinches off so quickly that it cannot be tra- 
versed even by a test signal moving with the velocity of 
light [isj . In order to prevent the shrinking of a worm- 
hole and to make it traversable, it is necessary to thread 
its throat with so-called "exotic matter", which is mat- 
ter that violates the averaged null energy condition (see 
@) Hj Many questions associated with the exis- 

tence of wormholes and with exotic matter remain un- 
solved, e.g. its stability against various processes, and 
different views have been expressed in the literature (see 
for example 0, [H, [H, [13 )• We will not discuss this 
here. Different types of wormholes may exist depend- 
ing on the type of exotic matter in their throats (see 
[i, [H, [T^). For example, it could be a "magnetic ex- 
otic matter" in which the main component is a strong, 
ordered, magnetic field plus a small amount of a "true 



exotic matter" (see [9| ) . Another type of exotic matter 
is a "scalar exotic matter" in the form of a scalar field 
with a negative energy density (see [3). Finally, it can 
be in the form of a mixed "magnetic-negative dust ex- 
otic matter" type in which it is a mixture of an ordered 
magnetic field and dust (matter with zero pressure) with 
negative matter density (see (l^)- The physical proper- 
ties of different types of wormholes are different (see for 
example 119]). The possible consequences of the passage 
of radiation through a wormhole with a "scalar exotic 
matter" was considered, using an analytical approach, 
by Doroshkevich, Kardashev, D. Novikov and I. Novikov 
in where a toy model was used to accept some very 
artificial hypotheses about the final state of the worm- 
hole. The results obtained in [131 allowed the authors to 
reach some important conclusions, but of course without 
numerical computations, the real dynamics of the pro- 
cesses of radiation propagation through a wormhole could 
not be considered. Other important analytical analyses 
of the existence and evolution of the wormholes and their 
possible transformation into black holes (BH) was seen 
in [2ll,[23,[2l[23,[2i,[2i,[li[28'|. Other important results 
are also seen in the papers [29, 30, .31). 

The new period of the investigation of the wormholes 
began from the seminal papers [33 . [3^ , where the worm- 
hole dynamics were analyzed using numerical simula- 
tions. Paper [33| is devoted mainly to the evolution of 
an initial static wormhole maintained by the scalar exotic 
matter under the gravitational action of the scalar radia- 
tion pulse passing the wormhole from one asymptotically 
flat region to other one. In paper [s^ the authors ana- 
lyze the nonlinear evolution of the wormhole perturbed 
by the scalar field. 

The goal of this paper is to continue the line of the 
investigation of these two works. We focus on the dy- 
namics of the scalar fields during the nonlinear stages of 
the wormhole evolution under the gravitational action of 
the ingoing pulse of the scalar radiation (both exotic and 
usual). These dynamics can explain some results of the 
papers [H, [11]. 
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This paper is organized as follows: 

In section |lT] we present the initial model of the worm- 
hole. In section IIIII the field equations are written out. 
The initial value problem is discussed in section IIVI In 
section |V] we present the numerical solution of the field 
equation which describe the passage of the signal through 
the wormhole and its subsequent evolution with the for- 
mation of a BH. In section IVII we analyze the case of 
a strong signal. The passage of the signal of the exotic 
scalar field is considered in section IVIII Finally we sum- 
marize our conclusions in section IVIIIl Mathematical 
details and the details of the numerical code are given in 
appendices IXllCl 



II. THE MODEL 

We study, numerically, the evolution of a simple 
spherical model of a wormhole maintained by an exotic 
scalar field and perturbed by an in-falling massless, self- 
gravitating scalar field. This initially in-falling scalar 
field is a pulse in the form of a spherical layer with some 
width and amplitude. It imitates a pulse of radiation 
directed into the wormhole and propagating with veloc- 
ity of light. We do not impose that the in- falling scalar 
field is weak. We will consider two cases: when this in- 
falling field has a positive energy density (imitating real 
radiation), and when it has a negative energy density 
(imitating exotic radiation). As an initial static worm- 
hole, let us consider the special case of the Armendariz- 
Picon solution of the Einstein equations for the spherical 
static state with the effective mass equal to zero [l^ ■ 
Another name of this wormhole is Moris-Thorne's (MT) 
wormhole. The corresponding line element for such a 
wormhole can be written as follows: 



space, but the wormhole (two asymptotically flat three- 
dimensional regions joined by a three-dimensional tun- 
nel) nevertheless exists due to the distribution ([2]) of the 
exotic scalar field. Of course there are not any apparent 
and event horizons in such an object and test signals can 
pass through the timnel in both directions. 

The initial motivations of our choice of the MT worm- 
hole as an initial static wormhole were the following. 

(1) The corresponding solutions (HJ-lH]) of the Einstein 
equations are very simple and elegant. 

(2) It was declared in [3| that the MT wormhole is 
stable. 



However the analytical |34l | and numerical j33i | anal- 
ysis (see also Appendix]^ demonstrates that the solu- 
tions ©-(111) are unstable against small spherical pertur- 
bations. Our numerical experiments confirms this con- 
clusion. 

In the light of this fact, the analysis of the possibility 
of the passage of the signals through the MT wormhole 
and their distortions have a special interest. 



III. FIELD EQUATIONS 

For the numerical analysis we will use the double null 
coordinates. The general line element in these coordi- 
nates can be written as: 



(3) 



where a{u^v) and r{u,v) are functions of the null coor- 
dinates u and V (in- and out-going respectively). The 
non-zero components of the Einstein tensor are: 
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= Q^ + R\ (1) 



where R is running from — oo to -l-oo, Q is a charac- 
teristic of the strength of the exotic supporting field ^ 
and dfi^ = dO'^ + siv? 9 dip"^ is the line element on the 
unit two-sphere. The general equations will be given be- 
low. Here we give the non-zero components of the energy- 
momentum tensor of the ^t-field for the solution ([T]): 
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The narrowest part of the wormhole (H]) is at i? = 0. This 
throat corresponds to r = Q- The physical analysis of the 
metric ([T]) is given in Here we emphasize that for 

the metric ([1]), in this reference frame, there is no gravita- 
tional acceleration at any point in the three-dimensional 
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The energy-momentum tensor can be written as a sum 
of contributions from the exotic scalar field and the 
ordinary scalar field $ with the positive energy density: 
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The u — u. V — V, u — V and 9 — 9 components of the 
Emstein equations respectively are (with c = 1, G = 1): 



r,uu - 2 r „ „ - r (*,„)' + r ($,„)' = (7) 
r „^ - 2 r „ cr,„ - r (*,„)^ + r ($,^)2 = (8) 



r,uv ^ 1- -z- = 

r Zr 

o2(T 



(9) 



r2 2r^ 



'f.u'f.v + = (10) 



The scalar fields satisfy the Gordon-Klein equation 
VV^* = and V''V^$ = 0, which in the metric ^ 
become: 

+ - (r^^-.n + r.„^',„) = (11) 
r 

$,™ + -(^i,$,« + r„$,,„) = (12) 



Equations (|9|l- (fT2l) are evolution equations which are 
supplemented by the two constraint equations ([7]) and 
([5]). It is noted that none of these equations depends 
directly on the scalar fields 5* and $, but only on their 
derivatives, i.e. the derivative of the scalar field is a phys- 
ical quantity, while the absolute value of the scalar field 
itself is not. Specifically we note the = -(*^„)2/(47r) 
and 7"* = — (5'^i,)^/(47r) components of the energy- 
momentum tensor ® and T*„ = {^^uY/{An) and T*„ = 
i^.v)^ / (47r) components of the energy-momentum tensor 
([S]), which are part of the constraint equations. Physi- 
cally Tuu and Tyy represents the flux of the scalar field 
through a surface of constant v and u respectively. These 
fluxes will play an important role in our interpretation of 
the numerical results in sections |Vl IVII and IVIII 



IV. INITIAL VALUE PROBLEM 

We wish to numerically evolve the unknown func- 
tions r{u,v), a(u,v), <I>(m, u) and ^'(u,w) throughout 
some computational domain. We do this by follow- 
ing the approach of [H [H, 113, H [H, SU, IH, 
|43| to numerically integrate the four evolution equa- 
tions ([5|)- P^ . These equations form a well-posed ini- 
tial value problem in which we can specify initial val- 
ues of the unknowns on two initial null segments, 
namely an ingoing [v — va — constant) and an outgoing 
{u = uq = constant) segment. We impose the constraint 
equations ([7]) and ^ on the initial segments. Consis- 
tency of the evolving fields with the constraint equations 
is then ensured via the contracted Bianchi identities [s^] , 
but we use the constraint equations throughout the do- 
main of integration to check the accuracy of the numeri- 
cal simulation. 

Our choice of the initial values corresponds to the fol- 
lowing physical situation; There is an MT-wormholc and 
at some distance from one of the entrances, at the initial 



moment there is a rather narrow spherical layer of an 
in-falling scalar field. There is not any radiation coming 
to the wormhole from the side of the other entrance. As 
mentioned, we specify the initial values on two initial seg- 
ments, namely u = constant and v ~ constant. In the 
next section |V] we consider the case where the in- falling 
layer consists of the <I>-field. Later in the section IVIII we 
consider the case of the in-falling layer consisting of the 
^--field. 

For the case of the layer from the $-field, the initial 
condition could be specified as follows: 

First, let us consider the static MT solution 



{MT) 



r2 Q2 + ^2 



(MT) _ 



arctan 



(13) 

In u — V coordinates we can write the MT solution as: 



r{u, v) = \j Q"^ + '^{v ~ u)^ 



^(A/T) ^ arctan 
$(^^^) = 



2Q 



(14) 
(15) 
(16) 
(17) 



Equations - (fTB|) completely specifics the MT- 
wormhole inu — v coordinates in all of our computational 
domain including the initial surfaces u — constant and 
V — constant. 

Now, let us consider the more general case with non- 
trivial <& and 5" scalar fields. In sections |V] and IVII we 
consider the case of a non-zero 4> field and in section 
IVIII we consider the case of a perturbed MT 5" field on 
the initial outgoing u = constant surface. In general we 
are free to choose the $ and ^ fields (and their deriva- 
tives) on the initial surfaces in any way we wish. Our 
choices for the $ and ^f-ficlds are described in sections 
IVl- IVIII We are also free to choose r{u,v) on the initial 
surfaces, this merely expresses the gauge freedom asso- 
ciated with the transformation u — > u{u),v v{v) (the 
line element ([3]) and the equations ([7)l- P^ are invariant 
to such a transformation). We continue to use eq. (fH|) 
on the initial surface as our choice of gauge. Hence, the 
only variable left for us to specify on the initial surfaces 
is cr. This can easily be found by integrating the con- 
straint equations eq. ([7|) and ([5]), which ensures that the 
constraint equations are satisfied on the initial hypersur- 
faces. Specifically on the outgoing u — uq hypersurface 
it is found by the integral: 



r (^^.)'+r (^,.)^ 
2r„ 



dv 



a{uo,v) = In ^-i= 

u 

(18) 

The initial conditions on the initial ingoing hypersur- 
face is set equal to the MT solution for all simulations. 
Hence, by specifying a distribution of the scalar fields <& 
and \1/ on the initial null segments, choosing a gauge and 
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charge parameter "Q" we can specify complete initial 
conditions on the initial null segments. Using a numeri- 
cal code (described in appendix |B| , we can then use the 
evolution equations, eqs. ((9|)- (fT2|) to evolve the unknown 
functions throughout the computational domain. 



V. PHYSICS OF THE PASSAGE OF THE 
'I'-FIELD PULSE THROUGH THE MT- 
WORMHOLE 

We investigate the full nonlinear processes arising in 
the case of propagation of a compact pulse of the scalar 
field through an MT-wormhole, using results of our nu- 
merical simulations. Our numerical code is described in 
appendix [B] and tested in appendix [Cl In this section we 
analyze the passage of the $-field pulse. In section IVTIl we 
investigate the passage of the ^J^-field pulse. Throughout 
this paper, the constant Q (see eq. ([Ij), prior to influence 
from scalar pulses, has initial value of Q = 1. The com- 
putational domain for all results throughout this paper is 
V = [8, 28] and u = [0, 20]. The flux of the scalar $-field 
into the wormhole is specified along initial u — = 
outside the wormhole in the following way: 



$ „(-uo,i;) = v4* sin^ 



Vl 



V2 - Vl 



(19) 



where vi and V2 marks the beginning and end of the ingo- 
ing scalar pulse, respectively, and A* measures the am- 
plitude of the pulse. Before and after the pulse, at w < wi 
and V > V2 the flux through u = uq is set equal to zero, 
i.e. ^^^(uo, v) = 0. The flux of the scalar <I>-field through 
initial segment v = vq is set equal to zero: $^„(u, vq) — 0. 
The expression p9)) can readily be integrated to give: 



A* 



27r {v — Vl) — {v2 — Vl) I 



27r 



V — Vl 



V2 - Vl 

Note that we formulate the initial condition directly 
for the fiux Ty^ = ($^„)^/(47r) of the scalar field through 
the surface u = uo, rather than for <I> itself since Tyy has 
the direct physical meaning. Also remember from section 
llVl that once the flux through the two initial surfaces has 
been chosen, all other initial conditions are determined 
by the model. In the examples of the results of our com- 
putations we specify vi = 9 and V2 = 11- In our compu- 
tations we vary the amplitude A* of the pulse. 

In Fig. [1] is seen the static case of the MT metric ([1]) 
in u — V coordinates, without any perturbations of the 
scalar fields. The throat r = 1 corresponds to the diago- 
nal u = V. 

In Fig. [5] one can see a typical example of the evolu- 
tion of the wormhole due to the action of the in-falling 
$ scalar field. Initial values for this <I>-field are vi = 9, 
V2 — 11, A* = 0.01. Fig. 2(a) shows the evolution of the 
lines r — const. As is seen in Fig. 2(a)[ the evolution cor- 
responds to contraction of the throat down to r = 0, and 
contraction of the whole wormhole. The real singularity 
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FIG. 1: Static MT-wormhole. Plot shows contour lines of 
constant r. Thick diagonal line "a" marks the throat of the 
wormhole (corresponding to the diagonal u = v) with r = 1. 
The thin lines marks lines of constant r increasing outwards 
from the throat of the wormhole with a constant spacing be- 
tween lines of Ar = 0.5 (i.e. line "b" marks r — 1.5 up to 
line "c" at r = 12.). 



of the space-time arises at r = 0. It is seen from the fact 
that when we come to r = 0, the metric coefficient e^*^ 
in formula ([3]) tends to infinity (see Fig. |2(b)"| ). Also ap- 
parent and event horizons arise. The apparent horizon 
corresponds to events where world lines r = const are 
horizontal or vertical (see 4j| page 85). All space-time 
is divided into R- and T-regions (see [l^, IH, S, • 
the _R- regions (lines "a" and "b"), the lines r = const are 
time-like and signals can go both to higher and smaller r. 
In the T-region (between the two apparent horizons) the 
lines r = const are space-like and a signal can propagate 
to smaller r only. The structure of the apparent horizon 
depends on the value of the fiux of the energy through 
it. We will analyze it later (see section IVTl) . 

Now we consider the evolution of the <&- and ^'-fields. 



Fig. 3(a) represents the evolution of the flux of the 
scalar $-energy into the wormhole. Fig. |3(b)| shows the 
evolution of the T*„ flux which arises as a result of T*, 
being scattered by the space-time curvature. Fig. 4(a 
represents the evolution of the difference of the — 
of the fluxes of the ^f-fleld in and out of the wormhole. 
The reason why for the ^f-fleld we show the difference of 
the in and out fluxes rather than the and T^-fluxes 
separately, is the fact that even for the static MT-solution 
(P) , (21) , the values of and are not equal to zero 
(but equal to each other). So only the difference between 
them is important for the dynamics. Finally Fig. |4(b)| - 
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(a) Lines of constant r, ranging from r = to r = 5 with a spacing 

between lines of Ar = 0.25 (line "a" marks r = 1.25, line "b" 
marks r = 5 and line "c" marks r = 0.75). The thick dotted lines 
marks the position of the apparent horizons. The thick fully drawn 
line marks the position of the singularity r = 0, hence the region 
marked "X" is outside of the computational domain. 



(b)Plot of the metric coefficient e^"^'^' as a function of distance 
parameter L = \/ (v — vq)^ + (u — uq)^ measuring the distance 
along the wormholc throat at it = u (with L = corresponding to 
where the <E>-pulsc first reaches the wormhole throat, i.e. 
vo = uo = 9). 



FIG. 2: Results of simulation with a non-zero <l>-field modeled after eq. (|20p with ui = 9, ?;2 = H and amplitude A"^ 



0.01. 



|4(d)| shows the comparison of all three fluxes at different 



slices: 



3, 14, 20. 



In Fig. 



3(a) and 3(b) the initial pulse (between 
9.0 < V < 11.0) and subsequent tails with resonances are 
clearly seen. In different regions, and are con- 
verted into one another due to curvature scattering and 



resonances. But at later u (see Fig. 4(c)[ |4(d)[ ) and for 
w > 11 the flux T:* - r,* of the ^--field dominates abso- 
lutely (in modulus). Thus for these regions this flux of 
the ^'-ficld determines the dynamics of the wormhole. 

From the Fig. [2]|4] the general picture of the evolution 
of the wormhole under the action of the passage of the 
compact $-field signal looks like the following. At the 
beginning, a rather weak signal produces small perturba- 
tions of the static wormhole. These perturbations trigger 
the evolution of the ^-field that maintains the wormhole. 
The fluxes of the ^-field (in- and out-fluxes) determine 
the subsequent evolution of the wormhole and lead to 
its collapse. At the beginning of this process the fluxes 
of the ^f-field are directed to the throat from both sides 
of the wormhole (regions A and B in fig. 4(a)). This 



flux to the throat of the ^'-field with the negative en- 
ergy density may lead to the formation of the negative 
local effective mass. Remember that before the begin- 
ning of the process the local effective mass was equal to 



zero everywhere (no any gravitational forces). Just af- 
ter the passage of the perturbing compact ingoing signal 
of the <I>-field, the effective mass become positive by the 
mass of the ingoing <I>-shell. In the process of the col- 
lapse the fluxes directed in the opposite directions (out 
of the throat) arise (regions C and D in fig. 4(a) ). These 



fluxes partly propagate out of the openings of the worm- 
hole and partly propagate inside of the forming BH. The 
fluxes of the 'I'-field going out of the wormhole carry away 
negative energy providing origin of the positive mass of 
the forming BH. That is emphasized in papers (3^. [33|. 
Everything looks like a collapse and an explosion of the 
*-field. 



A very important fact is that the compact signal of 
the $-field propagates through the wormhole from one 
asymptotically flat region into another asymptotically 
flat region before the beginning of the collapse, this is 
especially clear from fig. |3(a) Only a small part of the 
<I>-signal is scattered back and into the arising BH. After 
the BH formation the propagation of a signal from one 
asymptotically flat space into another one is impossible 
in any direction. 
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(a) Contour lines of the ingoing $ flux, = const.. The hnes "a", (b) Contour hnes of the outgoing $ flux, r*„ = const.. The Unes 
"b" and "c" denotes regions of = 10~^,T^^ = 10~^ and "a", "b" and "c" denotes regions of T*„ = 3 ■ 10~^,T*^, = 3 • 10~^ 
= 3 ■ 10~^ respectively. and T*„ = 10~® respectively. The letter "d" identifles a region 

where the outgoing flux increases towards the r = singularity, the 
contour line closest to the r = singularity corresponds to 
T* = 10"" 

FIG. 3: Contour lines of the in- and out-going <E> fluxes for the simulation with v\ = %,V2 = 11 and amplitude A* — 0.01. In 
both plots lines are logaritmically spaced with a factor 10 2 between the lines. The thick solid line marks the location of the 
r = singularity and the region "X" marks the area outside of the computational domain. 



VI. STRONGER $-SIGNAL, HORIZONS AND 
RATE OF THE COLLAPSE 

Now let us consider the case of an esse ntially stro nger 
infalling $-signal with A* = 0.5 (see Fig. |5(a) 5(d) ). 

In this case, the collapse of the wormhole arises much 
faster. The essential part of the initial (f>-signal comes 



into the arising BH, this can be seen from figure 5(b) 
The picture of the process is very asymmetrical. 

The asymmetry is clearly seen at the picture of the 
apparent horizon (figure 5(a)). In general, with any 
strength of the compact $-signal, the apparent horizon 
arises at the throat just at the moment of the passage 
of the <I>-signal (see Fig. 2(a) ). At the beginning the dy- 



namics of the wormhole is rather slow (if the $-signal is 
weak enough) and two branches of the apparent horizon 
are very close from each other. 

When the fast collapse begins these branches are going 
in different directions. 

In the case of a strong initial pulse of the <I>-field 
(Fig. 5(a)) the shape of the branch of the apparent hori- 
zon from the side of the coming <I>-field signal is quite 
different from the one at the opposite side. 

The analogous change of the apparent horizon un- 
der the action of the incoming signal we observe in the 



charged BH, see |42| . 

The picture of the evolution of the ^'-field in the case 
of the strong signal (see Fig. 5(d) ) is qualitatively similar 
to the case of weaker signals (Fig. 4(a) I. 



About the event horizon. From Fig. 2(a) it is clear that 



events close enough to the singularity r = are inside of 
the BH and thus inside of the event horizon. Indeed, 
the light signal from these events (vertical and horizon- 
tal lines) will come to the singularity and cannot go to 
infinity. At Fig. 2(a) event horizon practically coincide 
with vertical and horizontal asymptotes to the singular- 
ity r = 0. The analogous region inside the BH can be 



found at the Fig. 5(a) 



Now let us consider the process of the collapse of the 



throat. In Fig. 6(a) one can see this process for different 
amplitudes A* of the ingoing signal of the $-field. It is 
seen that for small A"^ the wormhole practically does not 
responds to the <I>-signal during long periods and after 
that collapses very fast. It looks like a very nonlinear 
process. 



Fig. 6(b) represents the dependence of the moment of 
collapse of the throat on the amplitude A* of the ingoing 
^-signal. It is seen that for small A* the period between 
the passage of the signal and the moment of the collapse 
can be very long. 
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(a)Contour lines of the resulting 'I'-fiux = — with a factor (b)Comparison of the <I> and ^ fluxes along the line u = 3. Line 



of 10 2 between lines. The lines "e" and "f denotes 
TZ - = +10-'^ and T„* - T*„ = -10"'' respectively, with the 

resulting absolute flux increasing towards the r = singularity. 
The regions marked "A" and "C" marks — > 0, "B" and 
"D" marks regions with r„* - T*„ < 0. The thick dotted line 
marks the contour line — = and the fully drawn thick 
line marks the r = singularity as in previous plots as well as the 
region "X" which marks the area outside of the computational 
domain. 



"A" marks | — T^^ | , line " B" marks T!^^ and line " C" marks 
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(c) Comparison of the 3> and VP fluxes along the line u = 14. Line 
"A" marks |X„* -T*„|, line "B" marks r*„ and line "C" marks 



(d)Comparison of the $ and VP fluxes along the line u = 20. Line 
"A" marks |r*„ - T*„|, line "B" marks r*„ and line "C" marks 



FIG. 4: Results of the simulation with vi =9,V2 = 11 and amplitude A* = 0.01. 
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(a)Lines of constant r, from r = to r = 5 with Ar = 0.2 (lines (b)Lines of constant T^. Lines are logarithmically spaced with a 




10 12 14 16 18 20 8 10 12 14 16 18 20 



(c)Lines of constant T^. Lines are logarithmically spaced with a (d)Lines of constant — T^, Region "A" and "C" marks 
factor 105 between lines, line "a" marks T* = lO-^-S. regions with T*„ - T*„ > 0, "B" and "D" marks regions with 

Tvv ~T^u<0 and the thick dotted line marks T„* - r*„ = 0. 



FIG. 5: Results for a simulation with vi = 9,V2 = 11 and amplitude A* = 0.5. For all figures, the thick fully drawn line marks 
the position of the singularity r = 0, and the region marked "X" marks area outside of the computational domain. 
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(a)r as a function of the distance parameter 
L = \/ {v — vo)^ + (u — no)^ with ijq = tto = 9- The curves "a" 

"b", "c", "d", "e", "f and "g" marks simulations with 
A* = 0.5,0.2,0.1,0.01,0.001,0.00035 and 0.0001 respectively. 



(b)Distance parameter L = \/{v — vq)'^ + (u — uq)'^ with 
vo = uo = 9 measured along u = v from the start of the pulse to 
the position of the r = singularity as a function of scalar field 
amplitude A*. 



FIG. 6: Development of the r — Q singularity as a function of pulse amplitude (with identical pulse widths of vi = 9, U2 = 11). 



VII. PASSAGE OF THE *-FIELD 

In this section we briefly consider the passage of the ^'- 
field compact pulse through a MT-wormhole. As we told 
at the end of section [ill the MT wormhole is unstable 
against small spherical perturbations. Here we consider 
not small perturbations in the form of a compact signal 
coming from outside. As we will see, the evolution is in 
agreement with our conclusion about instability of the 
MT wormhole. Now we put the <I>-field= everywhere 
and model the ingoing ^'-flux as a perturbation to the 
background vpfA^"^) field supporting the wormhole: 



(21) 



The form for the perturbation of the ^f-field into the 
wormhole is specified almost the same way as in the case 
of the $-field in the section |Vl But now instead of the 
formula (1191) we should write down 



■0,1,(^0, v) = A sin TT 



V — Vl 
V2 — Vl 



(22) 



All other remarks about the initial conditions from sec- 
tion |V] are correct here also. But now there is one new 
factor in this case. Namely, this pulse of the ^P-field is an 
addition to the background of the Vf-fleld of the MT so- 
lution. This addition can be taken with the sign plus or 
minus. This means that A* can be positive or negative. 
The evolution depends critically on this sign. 

Fig. [7] represent the case of the expression (|22l) with the 
sign A* > 0. The wormhole expands (Fig. 7(a) I and the 



size of the throat becomes larger and larger. The physical 
reason for this expansion is clear, the additional amount 
of the \E'-field with the negative energy density creates 
gravitational repulsion. This gives the initial push. The 
expansion at the late w-parameter looks like an almost 
"inertial". 

Between two branches of the apparent horizon any test 
particle or radiation can move to greater and greater r 
onl y. T his is a so called expanding r+ region, see 

Fig. |7(b)] represent the fluxes of the ^-field out of the 
throat to greater r. 

The case of the expression ([2^ with the sign < is 
represented in Fig. [51 Now there is a deficit of the ^f-field 
as compared with the static solution ((1]), It leads to 
collapse of the wormhole. Qualitatively the process of 
the collapse is the same as it was described in section [Vl 



VIII. CONCLUSIONS 

In this paper we investigated numerically the process 
of the passage of the radiation pulse through a worm- 
hole and the subsequent evolution of the wormhole as 
respond to a gravitational action of this signal. It is the 
continuation of the investigations started in the papers 
[s^ . [ssl . [b^ I . In our work we focus mainly on the evo- 
lution of the fields $ and 5". We use the MT solution 
[see ((I])-(ll|)] of the Einstein equations as a model of an 
initial static wormhole. 
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(a)Lines of constant r, from r = 1 to r = 5 with Ar = 0.5 between (b)Lincs of constant — T^^, Region "a" marks region with 
lines (Une "a" marks r = 1.5 and Une "b" marks r = 11, — < "b" marks region with — > and the thick 

dotted line marks T* — T* = 0. 



FIG. 7: Results for a simulation with vi =9,V2 = 11 and >1* = +0.01 ($ = everywhere). 




(a)Lincs of constant r, from r = 1 to r = 10 with spacing Ar = 0.5 (b)Lincs of constant — T^, Regions "A" and "C" marks 
between lines, lines "a" marks r = 1.5, line "b" marks r = 10. regions with T„* - T*„ > 0, "B" and "D" marks regions with 
Thick dotted lines marks positions of apparent horizons. — < and the thick dotted line marks — = 0. 

FIG. 8: Results for a simulation with vi = 9, 02 = 11 and A* = —0.01 ("1> = everywhere). As previously, the thick fully 
drawn line marks the position of r = and the region marked "X" indicates the region beyond our computational domain. 
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The spherically symmetrical radiation pulse was mod- 
eled by the self-gravitating, minimally coupled, massless 
scalar fields $ and ^ [see eqs ([T9l) -(|22 |) ]. We created and 
tested a numerical code which is stable and second or- 
der accurate. For our computations we used an adaptive 
mesh refinement approach in both ingoing u and outgoing 
V directions. We have demonstrated that if an amplitude 
of the $-field, A*, is small enough, then just after the 
propagation of it through the throat the perturbations 
of the static wormhole are small. The compact signal of 
the $-field propagates through the wormhole from one 
asymptotically flat region into another one. The scatter- 
ing of the signal is small also. But the small perturbations 
trigger the evolution of the ^'-field which maintains the 
wormhole. The fluxes of the ^f-field (in- and out-fluxes) 
determine the subsequent evolution of the wormhole and 
lead to its collapse and a BH arises. 

After the BH formation the propagation of a signal 
from one asymptotically flat space into another one is 
impossible in any direction. 

In the case of the strong enough <i>-signal in-falling into 
the wormhole, the collapse of the wormhole arises much 
faster. The essential part of the initial <i>-signal falls into 
an arising BH. 

We analyzed the dependence of the period between the 
passage of the signal through the throat and the moment 
of the collapse of the throat in r = on the amplitude 
A*. For smaller and smaller amplitude A*, this period 
becomes longer and longer. 

Note, that the collapse and the BH formation after the 
passage of the $-signal through the wormhole is typical 
for this special type of the wormholes supported by a 
'3/-field. Other types of wormholes may behave quite dif- 
ferent. We will investigate such cases in another paper. 

Finally we investigated the passage of a compact signal 
of the \E'-field through the wormhole and evolution of the 

field in this case. 

As we emphasized above our work continues the work 
[H,!!^! . We investigate mainly the evolution of the $ and 
'3/ fields during the nonlinear evolution of the wormhole. 
On the other hand we want to emphasize the difference in 
our approach from the approach in the paper [33| . In [33| 
the authors formulate the initial condition for perturba- 
tions in such way that it corresponds to two fluxes of the 
scalar fields which propagate in positive and negative di- 
rections. The center of the perturbations is at the throat 
of the wormhole. Our formulation of the initial condition 
corresponds to a compact signal coming into the worm- 
hole from the right-hand asymptotically flat region. This 
relates our approach with possible applications to worm- 
hole astrophysics. 



our attention to their important papers. 

This work was supported in part by the JSPS 
Postdoctoral Fellowship For Foreign Researchers, the 
Grant-in-Aid for Scientific Research Fund of the JSPS 
(19 — 07795), Russian Foundation for Basic Research 
(project codes: 07 - 02 - 01128 - a, 08 - 02 - 00090 - a, 
08 - 02 - 00159 - a) , scientific schools: NSh - 626.2008.2, 
Sh — 2469.2008.2 and by the program Origin and Evolu- 
tion of Stars and Galaxies 2008 of Russian Academy of 
Sciences. 

The authors thank the Niels Bohr Institute for hospi- 
tality during their visit. 



APPENDIX A: INSTABILITY 

For the spherical metric 

ds^ = -dr^ + e^dR^ + d^^ , = e''(i?^ + Q^).(Al) 



the Einstein's tensor is: 



GZ=e' 



(A2) 



i?(377-A),„ 2Q2 + i?2 



Y2 



y4 



y2 4 



I I RV,R 



(A3) 



Y2 '^'-^-^ 4^,T 



Go — 



y4 



(A4) 



mv-^),R 

2Y^ 



G" 



%rR + ^ U.H + ^ 



(A5) 



Here Y"^ = + Q"^ . 

The equation for the scalar field is 
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exp 1 '7 + 2 ) 
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exp [V^ I Y-'^ 



(A7) 
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The energy-momentum tensor is: 



(A8) 



rpti 



The static MT solution is 



Tj^O, A = 0, 



y2 



(A9) 



r],xx H ^ h 2- 



1 + a;2 1 + a;2 

and using one of the solution (jAlSp we get 



A = -277, ^,tt = V,xx 

For ?7 (X expiiijjt) we get 
6 a; 



rV.x 



1 + 



X'^ + 1 



x^ + 1 



(A19) 



ri = (A20) 



(AlO) 



A 

2 



/ 



(1+^2)3/2^ 



f.x 



2VT 



(A21) 



1. Perturbations 

Further on we will use dimensionless coordinates x = 
R/Q, t = ct/Q and will consider linearized Einstein's 
equations: G\. = SttT^. Perturbed components of the 
energy-momentum tensor are dimensionless. From (jA2f 
^8|) for 77 ^ 1, A <C 1 we get for perturbations: 

6Tl = ~5T^ = ST^ = |^[2V^,. - A*,.] (AH) 
an 



(l-ha;2)2 



/ = 



(A22) 



As is seen from equations (IA20p - (|A22l) functions r?(a;), 
f{x) and ip{x) include even and odd modes. This means 
that for the even mode the function 'ip{x) is extreme at 
x = while for the odd mode 'ip{0) = 0. 

The Eq. (|A22p is identical to the one-dimensional 
Schrodinger's equation with the energy E — llP' and the 
potential box 



5T 



(A12) 



U 



(l + a;2)2' 



Udx 



-l.Svr w -4.7 (A23) 



a;2 -t- 1 



(A13) 



qHgI 



(A14) 



ixrj^x + rj xX,x x'^ + 2 



l2e^_R, Xr^.x+V X^X 



x^ + 1 



(x2 + 1)2 



(A15) 



(A16) 



Tj.xx ~ A,tt - iitt x{\^x - '2r],x) A 



2 2(2:2 _^ 1) (2.2 1)2 

From the condition SGl — ~SG^ we get 



(A17) 



1 + a;2 l-f a;2 

and the condition SGl = '^^e leads to equation 

A,tt + Ti,tt = (A18) 



For such equations the general structure of solution is 
well known. For example, for lu^ > we have. 



f — p'f^^ ) (^1^^'^^ forx<C 1 ( \0A\ 

■' 1 C2e-*'^^-|-C3e'"^ forx>l ^ ' 



what are the incident and reflected waves for x ^ 1 and 
the traveled wave for a; <C — 1. 

For the potential (|A23[) there is at least one solution 



with 



-a , 3 > a > and for finite f{x) at x 



±00 we expect to get 

ci exp(aa;) 



for x ^ —I 



/ ~ e"' <( C2 sin(V3 - a^x + cj)) for - 1< x < 1(A25) 
C3exp(— ax) forx 3> 1 

with = const and a discrete set of a = a^, i = 1, 2, ... 
and ^/3>a^>0. 

To estimate the expected principal values a we can use 
solutions of the Schrodinger's equation with two similar 
potential boxes, namely. 



U = Ui = 



and 



U^U2^ -3ch" X, 



forx < —1 

for - 1 < X < 1 (A26) 
forx > 1 



-00 < a; < 00 . (A27) 
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For both potentials (K26\MXT7\i we have / = Udx = 
-6 what is close to (|A23p . For the potential (IA26P the 
principal value is a = ai w 1.4 and for the potential 
(|M7l) it is a = a2 ^ (VTS - l)/2 w 1.3. These results 
indicate that for Eq. (|A22p we can expect to get also the 
principal value a « 1.3—1.4. In the case the growing 
solution of Eq. (jA22p is located at < Xb ^ 0.6. 

This means that there is an instability in the central re- 
gion of the wormhole {\R\ < XbQ) while at \R\ > x^Q the 
corresponding function / is exponentially small. These 
results demonstrate the instability of the MT solution in 
respect to the small perturbations of the field ^E". 



APPENDIX B: DESCRIPTION OF THE CODE 

In this appendix we describe our numerical code used 
to obtain the results presented in the paper. 

To numerically integrate the evolution equations, eqs. 
([51)- (fT^ . throughout some computational domain, we 
have created a new numerical code. The code has adap- 
tive mesh refinement (AMR) capabilities in both in- and 
outgoing directions (contrary to our previous code, which 
had only AMR along the ingoing it-direction [i^l). The 
numerical scheme used to evolve the computational cells 
is identical to the one we used in [12] (which evolves the 
unknown variables with second order accuracy), while 
the AMR algorithm in our code is largely based upon 
the design of the algorithm presented in ^13|. 

A key feature of this AMR algorithm is the usage of 
a self-shadow hierarchy to compute the truncation error 
estimates which determines whether or not a computa- 
tional cell should be split to a higher level. A self-shadow 
hierarchy works by always simultaneously evolving two 
grids, one twice the resolution of the other. Wherever the 
computational points in two grids coincide, it is possible 
to estimate the local truncation error (TE) by calculat- 
ing the difference between the twos solutions. This TE 
estimate is then measured against some limit TEmax and 
if the local truncation error is larger then that allowed, 
the computational cell is split to a higher level. 

However, the point-wise computation of the TE as de- 
scribed above is, in general, not the optimal way of com- 
puting the TE as the solutions to the wave-like finite- 
difference equations is in general oscillatory in nature and 
will tend to go to zero at certain points within the compu- 
tational domain [43| . Therefore, in practice, when calcu- 
lating the TE we average the point-wise TE over several 
cells to the past of that point. Where the splitting struc- 
ture allow it, we use a total of nine points to the causal 
past of a point to calculate the TE which determines 
whether a cell should be split or not. Furthermore, we 
calculate the TE for all four dynamic variables, if any one 
of these become greater then our accepted limit TE^ax, 
the cell is split to a higher level. 

When splitting a cell to a finer level, we need to per- 
form an interpolation of the cells on the more coarse level 
in order to obtain initial data at the refinement bound- 




FIG. 9: Convergence factors for xi'^) (see text for definition) 
along lines of constant u for non-AMR simulations with vary- 
ing resolutions. Lines are for resolutions a) A'^ = 5, b) A'' = 10, 
c) TV = 20, d) N = 40 and e) iV = 80. 



aries for evolving the cells on the finer level (see [42] and 
[4^ 1 for details). If a sufficient number of adjacent points 
are available on the coarse parent level we use cubic La- 
grange interpolating polynomials, otherwise we use linear 
interpolation. 



APPENDIX C: TESTING OF THE CODE 



We have tested our code and found it to be stable and 
second-order accurate. As a first, most basic test, we 
have verified that the code does indeed converge to the 
static MT solution. 

In this section, we demonstrate that the code is also 
converging and second order accurate when including the 
effects of a non-trivial infalling $ field. The initial con- 
ditions for the tests in this section are similar to those 
described in section |Vl i.e. the initial conditions for the 

field is set to the exact MT solution, and the form 
of the infalling $ field is modeled after eq. p9|) with 
vi =9,V2 = 11 and A* = 0.0035. The strength of this 

pulse is such that the wormhole begins its contraction 
and almost, but not quite, forms a black hole within our 
computational domain. Our computational domain is, as 
for all simulations in this paper, in the range v — [8, 28] 
and u = JO, 20]. 
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(a)||Ci,u|| (see text) for the constraint equation (0. 



(b)||Cuu|| (see text) for the constraint equation 



FIG. 10: Average values of the constraint equations along lines of constant v for non-AMR simulations with varying resolutions. 
Lines correspond to resolutions of a) iV = 5, b) iV = 10, c) iV = 20, d) iV = 40 and e) N = 80. 



1. Non-AMR convergence tests 

First we test the basic convergence properties of the 
code by performing a series of simulations with varying 
base resolutions but with the AMR algorithm disabled. 

To estimate the convergence rate along a ray of u = 
constant we calculate : 



X{Xn) = - ^n, i _ T I 

where x]^ denotes the dynamic variable x at the i-th 
grid point (along the given u — constant line) of simula- 
tion with base resolution N (where N = 1/ Au ~ 1/ Au). 
The n points over which we sum, is taken along an out- 
going ray of u = constant. Naturally, in order for this 
expression to make sense, we are only using those n points 
which all three resolutions have in common. If the code 
is second-order convergent, we would expect to see a con- 
vergence factor of X = 4. 

We calculate expression (|C1I) along lines of many u — 
constant and for many base resolutions N, such that we 
get a picture of the convergence properties throughout 
all of our computational domain for varying resolutions. 
Figure [H] shows the convergence rates for the dynamical 
variable a for a large number of resolutions. The figure 
shows that for small u (below u = 10), the code is sec- 
ond order converging. The line e (corresponding to the 
base resolution N — 80) gives the impression of less than 
second order convergence at small u, however a closer in- 
vestigation of this behavior has revealed that this is due 



to the solution in this region has reached such a high 
degree of convergence, that numerical roundoff errors in 
the numerical representation affects the calculation of the 
convergence factor. 

For high u we see that especially for low resolutions, 
the convergence is far from second order. This is mainly 
due to the solution in this region, for low resolutions, is 
diverging too much from the physical solution. Such be- 
havior is not uncommon for non-linear systems. However, 
as is seen, for increasing resolutions, even for high m, the 
convergence factor of the solution becomes increasingly 
higher, indicating that the numerical solution is converg- 
ing with second order accuracy. That a fixed resolution 
can be more than adequate in one part of the computa- 
tional domain and insufficient in another part is a strong 
indication of the need for AMR in our code. The proper- 
ties of the code with AMR enables is investigated in the 
next subsection. 

The other three dynamical variables has convergence 
properties similar to those of a. 

It is not enough, however, to demonstrate that our 
code is converging, we also need to demonstrate that it 
is converging to a physical solution. We do this by cal- 
culating the average value 

1 " 

11^^11 = -E^V (C2) 

i 

of the two constraint equations along lines of constant 
u for different resolutions. Figure [TUl demonstrates that 
both of the constraint equations is converging towards 
zero for increasing resolution, thus demonstrating that 
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(a) Convergence factor ^(r). 



(b) Convergence factor for 5((t). 
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(c) Convergence factor for 5(0)- (d)Convergence factor for ^(tp)... 

FIG. 11: Convergence factors ^(x) along line of u = 20 for simulations with AMR enabled, see text for details. 



our numerical solution is indeed converging towards a 
physical solution for increased resolutions. 



2. AMR convergence tests 

An important component of our code is its AMR capa- 
bilities. In this subsection we demonstrate the behavior 
of the code when the AMR algorithm is enabled. 

For each simulation we must specify a resolution of the 
underlying base grid and we must also specify a trunca- 



tion error acceptance limit (TEmax) which determines 
whether or not a computation cell should be split. For 
all simulations we allow for a maximum of 8 splitting lev- 
els (including the basegrid and it's self-shadowing level, 
see Appendix [B] and 43]). 

To investigate the convergence properties for the code 
with AMR enabled we did three simulations with the 
same initial conditions as in the previous section, varying 
the base resolution and TEj^ax to mimic the doubling of 
resolution from one simulation to the next as follows; 
The lowest resolution simulation has a base resolution 
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(a) Average of Cw along lines of constant u for varying TEmax-- (b) Average of Cuu along lines of constant u for varying TEmax- 

FIG. 12: Average residual of the constraint equations (along lines of constant u) for fixed base resolution N = 5, AMR enabled 
and varying TEmax- Lines are for a) TEmax = le — 8, b) TEmax = le - 9, c) TEmax = le - 10, d) TEmax = le — 11 and e) 
TEmax = le - 12. 



of = 5 (corresponding to a total of 100 gridpoints) 
and a TEmax of rjv, the second simulation has a base 
resolution of A'^ = 10 and a TE^asx of = Tff and 
the last simulation has a base resolution of iV = 20 and a 
TEmasx of T47V = TAr/16. The varying of TEmax mimics 
the doubling of the resolution from one simulation to the 
next (under the assumption of second order convergence) 
and in principle the splitting level hierarchies should be 
similar for all three simulations. In practice, however, we 
do not obtain identical splitting hierarchies because the 
splitting structure is determined by the TE estimate and 
this will not scale exactly as the leading order part of 
the actual truncation error, which decreases by a factor 
of 4 each time the mesh spacing is halved for a second 
order accurate scheme [1^. Nevertheless, this test gives a 
good indication of the behavior of the code and if the test 
showed non-convergent results it would indicate problems 
with the implementation. 

In this section, to calculate the convergence factor, we 
use data along m = 20 (i.e. the upper border of our 
computational domain). As this is the outer border of 
our computational domain, convergence here is a good 
indication of the behavior throughout our domain. Fur- 
thermore, this border is closest to the r ~ singularity 
(that would be formed within our domain if the initial 
<l>-pulse had been slightly stronger), hence it is the most 
demanding place in our domain to test convergence. 

Figure [11] shows the convergence factors : 



1 \x' 



N 



(C3) 



where denotes the dynamic variable x at the j-th grid 
point (along u = 20) of simulation with base resolution 
N (where N = 1/Au — 1/Au). Naturally, in order for 
this expression to make sense, we are only using those n 
points which all three resolutions have in common. If the 
code is second-order convergent, we would expect to see 
a convergence factor of ^ = 4. 

In the figure, we see that for small v, in general, the 
factors are greater than 4, whereas for larger w, in gen- 
eral the convergence factors are smaller. It is clear that, 
especially for r and a we are seeing unrealistically high 
convergence factors for early v, however, as noted above, 
it would be surprising to recover perfect second order 
convergence in this test. The important point of fig. 1111 
however, is that all the dynamic variables are converging, 
also when the AMR algorithm is turned on, and that this 
convergence is in good agreement with that expected. 

Finally we show the results of a series of simulations 
with a fixed base resolution of = 5 but varying the 
TEmax limits. In fig. [T^lis shown the average of each of 
the constraint equations along lines of constant u, i.e. : 



\xte\ 



''TBI 



(C4) 



(where x is the constraint equation, i denotes the i-th 
out of n point along a line of constant u) . 

We see that both the constraint equations are converg- 
ing towards zero (although, for this test is makes no sense 
to talk about the "convergence rate"). It is also visible 
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FIG. 13: Demonstration of splitting level structure. Numbers 
indicate the number of split levels including the base shadow 
hierarchy (i.e. "2" is most basic structure possible in our 
configuration). 

from fig. 1121 that for increasing u, the constraint equa- 
tions are increasingly violated. This is only expected and 
is caused partly by truncation errors being accumulated 



throughout the computational domain and partly by the 
solution getting close to the region where a black hole is 
formed and hence also the region where strong gradients 
in the dynamic variables puts larger requirements on the 
numerical code. 

Another important point related with an AMR algo- 
rithm is how smooth the splitting levels are. It is well 
known that when a higher splitting level is introduced, 
this usually introduce high frequency noise into the solu- 
tion near the borders of higher resolution [43], which, in 
the worst case, may ultimately lead to the crash of the 
code. Because of this, it is highly desirable to have as few 
splitting levels as possible and to have as few splitting 
borders as possible. Figure [T2] illustrates the splitting 
level structure for a typical simulation with base resolu- 
tion and TEmax equal to the settings used in our actual 
simulations in sections IVllVIIl 

It is seen that the splitting structures are very smooth 
with no sudden or sporadic jumps, indicating that our 
implementation of the AMR algorithm is very successful. 
In [431 it is suggested to use a small numerical dissipation, 
but we have found that this is not necessary with our 
code, for the simulations presented in this paper. 

Finally it should be noted, that for all results presented 
in this paper, we have performed a great number of sim- 
ulations with varying base resolutions and TEmax- This 
we did to be certain that the results on which we have 
based our conclusions, had converged to such a degree 
that we could trust the results to be representative of 
the underlying physics. 
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